Influence of Heavy Quark Recombination on the Nucleon Strangeness Asymmetry 
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The nucleon strange and anti-strange distribution asymmetry is an important issue in the study 
of the nucleon structure. In this work, we show that the heavy quark recombination processes from a 
perturbative QCD picture can give a sizable influence on the measurement of the nucleon strangeness 
asymmetry from charged-current charm production processes, such as the CCFR and NuTeV dimuon 
measurements. When the influence of heavy quark recombination is considered, a positive effective 
^'S'hr should be added to the initially extracted strangeness asymmetry S~ = f dx[s{x) — s{x)]x, 
supporting the strangeness asymmetry S~ being positive, which is helpful to explain the NuTeV 
anomaly within the framework of the standard model. 

PACS numbers: 12.38.Bx, 13.15.+g, 13.87.Ce, 14.20. Dh 



The nucleon strange and anti-strange distributions are 
important quantities in the study of the nucleon struc- 
ture, and a clear knowledge of them helps for a better un- 
derstanding of some related phenomena in experiments. 
The nucleon strange quark-antiquark distribution asym- 
metry ispredicted naturally by some non-perturbative 
models and a positive strangeness asymmetry S~ = 
J dx[s{x) — s{x)]x has been shown 0, H, 0) B IE 
to be a promising mechanism to explain the NuTeV 
anomaly [§'] within the framework of the standard 
model. Perturbative QCD at three-loops can also gen- 
erate a strangeness asymmetry [l^], however, the ob- 
tained magnitude is one order smaller to be relevant to 
the NuTeV anomaly. 

In the measurement of the nucleon strangeness asym- 
metry, some valuable works have been done, though no 
conclusive result has been reached. Since strangeness 
asymmetry should be a very small quantity in inclu- 
sive deep inelastic scattering (DIS) cross sections, it 
is difficult to be extracted precisely. However, (anti- 
)neutrino induced charged current charm production pro- 
cesses are quite sensitive to the (anti-)strange distribu- 
tion, and thus can provide valuable information on the 
strangeness asymmetry. CCFR and NuTeV dimuon mea- 
surements are of such experiments. Although 
earlier analysis of dimuon data did not show support of 
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the nucleon strangeness asymmetry \lc 
cent next to leading order (NLO) analysis of the NuTeV 
data with improved method does show some evidence of 
the nucleon strangeness asymmetry [iHj S~ = 0.00196 ± 
0.00046(stat) ±0.00045(syst) ±0.00128(external). Mean- 
while, global analysis have indicated positive strangeness 
asymmetry iMli3,such as the most recent work of Lai 
et al. [i3| , who include both CCFR and NuTeV dimuon 
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data sets in their analysis, and produce the allowed range 
of -0.001 < S*- < 0.005 at 90% confidence level. These 
analysis suggest that S*^ is likely positive. 

In this work, we aim at checking the measurement of 
the nucleon strangeness asymmetry by including a per- 
turbatively calculable QCD effect. We find that the 
heavy quark recombination processes can produce a size- 
able influence on the measurements with charged-current 
charm production process such as the CCFR and NuTeV 
dimuon measurements. 

Heavy quark recombination [3, H) 01 combines a 
heavy quark, e.g., c quark, with a light antiquark q (or c 
with q) of relative small momentum in the hard scatter- 
ing, and the (eg) subsequently hadronizes into a D meson. 
Refs. [l^, [i^l employ simple perturbative QCD pictures 
and explain the charm photoproduction asymmetry and 
the leading particle effect [2l|, |2§| successfully. In the 
following, we show how the heavy quark recombination 
influences the measurement of the nucleon strangeness 
asymmetry. 

The CCFR and NuTeV dimuon measurements have 
provided important information on the strangeness de- 
grees of freedom in the parton structure of the nucleon. 
These measurements both rely on the (anti-) neutrino in- 
duced charged-current charm production processes, with 
the leading order (LO) subprocesses being i^^^ + s{d) 
fi~ +c and + s{d) /i+ + c. The produced c(c) quark 
hadronizes and then decays partially into fj,~^ {^^ ) to form 
a second /u. The oppositely signed dimuon events in ex- 
periment are then recorded for analysis of the nucleon 
strange distributions. 

The CCFR and NuTeV experiments use iron as their 
target, and for simplicity, we take it as an isoscalar target. 
The strange quark antiquark distribution asymmetry is 
directly related to the difference between neutrino and 
anti-neutrino induced dimuon differential cross section 
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FIG. 1; (a) and (b) are diagrams for cq recombination into 
a D meson in neutrino induced process ((Sjl; (c) and (d) are 
diagrams for cq recombination into a D meson in anti-neutrino 
induced process (|3]). Thick lines are heavy quarks, and shaded 
blobs are D or D mesons. 



at LO, which can be expressed as [22 



d^dy 



d^dy 
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xiam-mwcsi' + ^^[dAo+Mmvccii' } , (i) 



where ^ is the light-cone momentum fraction of the struck 
quark and is related to the Bjorken scaling variable x 
through ^ = x{l + ml/Q'^), S = 2ME„ and y = u/E^ 
with Ey and v being the incident energy and the energy 
transfer in the nucleon rest frame, = 1 + Q"^ /M"^ 
and fc=l~ m^/2MEy^, and Be is the branching ra- 



tio for 



~X. The valence contribution in Eq. ^ 



is suppressed relative to strange contribution from their 
relative coefficients |VcsP ^ 0.9 and \Vcd\'^ ~ 0.05, thus 
this cross section difference of Eq. ([1]) is sensitive to the 
strange distribution asymmetry S^{^) = S,[s{^) — s{^)]. 
The heavy quark recombination processes as 



+ q^H +s{d) + D{cq), 



Vf, + q^ + sid) + D{cq), 



(2) 



(3) 



(diagrams in FIG. 1) can also contribute to the dimuon 
final states through D{D) decays: D — > /i+X and D 

X . These processes are possible to have sizable effect 
in the extraction of the nucleon strangeness asymmetry. 

The processes of ([2]) and ([3|) contribute to the differ- 
ence between neutrino and anti-neutrino induced dimuon 
differential cross section at higher order, which can be ex- 



pressed as 
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where d(7]j(^cq) denotes the cross section for the subprocess 
(l2|), which is identical to the cross section of subprocess 
(l3|) from charge symmetry. i?£)(cg) is the branching ration 
for D{cq) — *■ fi'^X. q denotes a light quark flavor from 
the nucleon, which could be u or d, and the D{cq) meson 
could be either a scalar ^Sq state or a vector '^Si state. 

Such a contribution as Eq. Q serves as an additional 
part in the extracted strange distribution asymmetry, 
such as the NLO analysis of the NuTuV dimuon data [l5|, 
because the recombination processes as FIG. 1 are not in- 
cluded in the ananysis. Thus the realistic strange distri- 
bution asymmetry ^j^^ilO should be the analysed result 
'^aiiaiy(0 minus the contribution from heavy quark re- 
combination processes, which is a negative quantity from 
Eq. (HI) since the nucleon structure ensures 9(0;)— g(a;) < 
for q = u,d. From Eq. ^ and Eq. ([4]), one gets 
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with 
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where ^S^^{S^) > since to minus a negative quantity is 
equivalent to plus a positive quantity. Thus the realis- 
tic strangeness asymmetry S~^^y = J d^S^^^^{S,) should be 
larger than the experimentally extracted value according 
to the contribution from the unaccounted recombination 
processes. 

Now we proceed to estimate the size of <5S'jjpj(^). We 
follow the method in Ref. [3] to calculate the heavy 
quark recombination process. For the color singlet ^Sq 
D[cq) production, the following substitution is made in 
the parton amplitude: 

Vj{Pq)ui{pc) Xgj^mJ+{-^c - mc)j5- (7) 

Then set pq — XqPc in the amplitude and take the limit 
Xq — > 0. Thus the amplitude for color singlet ^Sq state 
D{cq) production is (diagrams in FIG. l(a)(b)): 

l&-KGF0lsmc5inf+ 
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where = u{k)j^{l — j5)u{p) is the lepton current. 

The 6ij in Eq. ([7]) is the color factor for color-singlet 
state, which is replaced by V^T^j for color-octet state 
together with the nonperturbative parameter /+ replaced 
by pi = and ps = (/+)^ will appear in cross 
sections to characterize the probability for a color-singlet 
and a color-octet ^So{cq) state to hadronize into a state 
including a ^S'o state D{cq) meson. The subprocess cross 
section for state D(cq) meson production thus can 
be expressed as 

dcfD(cq) = da[cq{^ Sf))i] ■ pi + da[cqC-So)a] ■ Ps- (9) 

The da[cq{^ So)s] can be calculated to be different from 
da[cq{^So)i] by a single color factor of 1/8. Thus, daD{cq) 
of Eq. ([9]) can be expressed as 

dcFD{cq) = da[cqCSQ)i] ■ Pcff [cg(^S'o) -> D{cq)], (10) 

with peff = Pi + Ps/S- 

For ^Si state production of vector meson D*{cq), sim- 
ilar substitution as Eq. ([7]) with the 75 replaced by / is 
made in the parton amplitude, where e is the polarization 
vector for the '^Si state. Similar expression as Eq. (fTO]) 
can be obtained for the subprocess cross section of vector 
D*{cq) meson production, 

daD'(cq) = da[cqCSi)i] ■ pMcqCSi) ^ D*icq)]. (11) 

Physically, there may be spin-flipped transitions such 
as cq(^So) — *■ D*{cq). While we have neglected such 
transitions, partly because the calculation of charm pho- 
toproduction [l^ and the leading particle effect [20[ have 
both set psf = 0, and partly because the inclusion of these 
transitions will not greatly affect our result, since both 
D{cq) and D*{cq) meson will decay similarly to 

The flavor-changing transitions, such as cu — > (cd), 
are also neglected as Refs. because these tran- 

sitions are relatively suppressed in the large Nc limit of 
QCD, and also because the inclusion of such transitions 
will not affect our result notedly. 

The number of free parameters can be greatly reduced 
from symmetries of the strong interaction. As discussed 
in Ref . [19J , heavy quark spin symmetry implies 

pMcqCSo) ^ D{cq)] = pMcqCSi) ^ D*{cq)], (12) 

and SU(3) light quark flavor symmetry indicates, for ex- 
ample. 



Thus, for isoscalar target, the SS^^iiO of Eq. (O can be 
expressed as 



Thus, only one parameter is left: 



(13) 



Psm = PoslcdCSo) D+] ^ p,s[cd{^Si) D*+] 

= p,s[cuCSo) -> D"] = pMcu^Si) -> D*"]. (14) 



dx[uy{x) + dy{x)] 

3< 



da[cqCSo)i], da[cq{'S,)i] ^ 
"1 I 7n "2 



d£.dy 



d£,dy 



■ Ps 



(15) 



where h = (5^,+ +Bdo)/2 and 62 = (Bd-^ +Bd*o)/2. 

The subprocess cross section can be calculated straight 
forward from the parton amplitudes (with the parame- 
ter /+ extracted out for psm)- Since the subprocess is a 
2 — !■ 3 process, there are five independent variables in the 
subprocess cross section. From the symmetry of the scat- 
tered /i~ around the incident direction, four independent 
variables are left, where two variables are transformed to 
^ and y (or Q^) and the other two are integrated out. 
Thus for 5S^Y{.i^) of fixed Q^, the integration is totaly 
of 3 dimensions including the integral on x in Eg. psp . 
The boundaries of the integration are determined from 
the allowed physical phase space. 

We use the CTEQ6L parton distributions for the nu- 
cleon [23 |. and the running coupling constant ag is as 
specified in CTEQ6L. We take TOc = 1-5 GeV and set 
the factorization scale to be -I- m^. Since the two 

muons in NuTeV experiment are required to have energy 
greater than 5 GeV, we try similar cuts for the produced 
p and the charmed meson in our integration. We find 
that the cross section from heavy quark recombination 
process decreases very slowly with the increase of the cut 
on the energy of the produced charmed meson. Thus the 
recombination processes are not suppressed by the cuts 
in experiments. 

FIG. 2 shows our result of 5S^^{£,) for = 160 
GeV, Q2 = 20 GeV and p^m 0.15. Such and 
are approximate averaged incident energy and in the 
NuTeV dimuon experiment [l3| . psm is the nonpertuba- 
tive parameter for the heavy quark recombination and 
Psm = 0.15 is the LO fitted result from charm photopro- 
duction asymmetry [l9l |. The branching ratios and \Vcs\ 
are taken to be the central values from Ref. [2^ . 

From FIG. 2, one sees that SS^^iiO is a valence- like 
distribution, with its peak in the range of ^ = 0.1 '--^ 0.3. 
From Eq. ([5|), the realistic strange distribution asymme- 
try is the sum of the analysed result of dimuon exper- 
iments and the effective contribution from heavy quark 
recombination <5'S'hr(C)- We can estimate 5S^^ by in- 
tegrating SS^Yiii) O'^sr ^, and get SS^^i ~ 0.0023 for 
Psni = 0.15. 

Such a value of 5S^^ significantly enhances the mea- 
sured strangeness asymmetry to a larger positive value, 
since S~^^y = S~^^^^ + SS^^i- Recent NLO analysis of the 
NuTeV dimuon data provides positive strangeness asym- 
metry centered at 0.00196 [15|. With the correction of 
the heavy quark recombination, the central value of the 



realistic strangeness asymmetry could be 



real 



0.0043. 
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FIG. 2: for = 160 GeV, = 20 GeV^ and 

Psm = 0.15. The dashed curve is the contribution from ^ So 
state; the dotted curve is the contribution from state; and 
the sohd curve is their sum, the iS^-^i^). 



Such a value of the strangeness asymmetry can explain 
the NuTeV anomaly to a large extent. NuTeV anomaly 
arises from the large discrepancy of the NuTeV measure- 
ment of the sin^u, with the standard model prediction, 
and becomes a hot debated area in recent years. The 
NuTeV measurement of the sin relies on the hypoth- 
esis that strange and anti-strange distributions are sym- 
metric. When this assumption is violated, their result on 
sin6'.uj will change. The influence of nonzero 5'~(^) to the 
result of sin dw is most sensitive in the range ^ — 0.06—0.3 
(FIG. 1 in Ref. [H), and such range is just the position 
of the peak of &S^^ (^) . Thus the non- vanishing 5S^^ (^) 
and S'~jjj may have a large effect in the NuTeV mea- 
surement of sin&uj. In such sensitive range, a positive 
strangeness asymmetry S~ of the order 0.005 can fill the 
gap between theory and the experiment of the NuTeV 
anomaly. The value of (55jjj^ w 0.0023 for ps,n = 0.15 
alone can provide nearly half of the strangeness asymme- 
try needed to explain the NuTeV anomaly. 

In our calculation, there are some uncertainties in the 
choice of the factorization scale [i and the parameter vfic- 
The impact for different choice of rric on our result is 
small (within 5 percent for mc = 1.5 ± 0.3 GeV). While 
for different choice of the factorization scale, the influence 
is large. Ref. [l3| calculate charm photoproduction with 
factorization scale yp^^+rof , where p_L is the transverse 
momentum of the produced D relative to the incident 
photon direction. In the process of this work, the photon 
is replaced by the W boson, and thus the calculation we 
take is performed for /.to = -|- m^, where is the 

transverse momentum of the produced D relative to the 
W boson direction in the nucleon rest frame. The result 
nearly trebles when /i = /io/2, and the result reduces 
nearly by half when /i — 2fiQ . The uncertainties may im- 
ply that higher order effects are still important. Ref. [l3| 
reports small effect on the predicted asymmetry when 
varying the factorization scale, where only the ratios of 



the cross sections are concerned, and scale dependence 
might be canceled in that case. 

The parameter psm still has some uncertainty in its 
value. In FIG. 2, we use psm = 0.15, which is from the 
fit of charm photoproduction asymmetry at LO [3] ■ As 
discussed in Ref. [19|, there is at least 30% uncertainty 
in this parameter psm due to finite heavy quark mass, 
SU(3) breaking and 1/Nc corrections, and more over, psm 
should be multiplied by a i^T factor if NLO corrections in 
photo-gluon fusion are incorporated in their calculation. 
Thus, parameter psm could well be as large as 0.3. Such 
a value for psm means SS^^ ~ 0.0046, which alone is 
sufficient to explain the NuTeV anomaly. On the other 
hand, psm could be smaller than 0.15, such as the LO fit 
from the leading particle effect pi = 0.06, where pi 
is the parameter for color singlet state, which is differ- 
ent from Psm but their size should be compatible. If we 
take Psm = 0.06, we get SS^^ ~ 0.0009. Such a SS^^i 
alone is too small to explain the NuTeV anomaly, how- 
ever, this (^iSj^j^ could still shift the dimuon result to a 
larger positive value, and make the allowed range of 
entirely positive. More precision determination of the 
parameters for the heavy quark recombination and a re- 
analysis of the dimuon events with consideration of the 
heavy quark recombination processes will be helpful to a 
better knowledge of the nucleon strangeness asymmetry. 

Our work implies the significance of using heavy quark 
recombination mechanism i.e., a perturba- 

tively calculable QCD effect, to reveal the strangeness 
asymmetry. Let us recall that a previous LO analysis Q 
of nucleon strangeness asymmetry by NuTeV collabora- 
tion reported a negative value S~ = —0.0027 ± 0.0013, 
whereas their new NLO analysis [l3| gave a positive value 
centered at 0.00196 as mentioned above. The difference 
between the two values is 0.0047, which is of the same 
order as the SS^^ estimated in this work. This again 
supports our work to take higher order effects into ac- 
count. 

In summary, we investigated the influence of heavy 
quark recombination in (anti-)neutrino induced charged 
current charm production processes on the measurement 
of the nucleon strange distribution asymmetry. Our re- 
sult shows that the influence could be quite sizable and 
the realistic strangeness asymmetry S~ should be larger 
than the initially experimental results. From our inves- 
tigation and the result of recent experimental analysis, 
the nucleon strangeness asymmetry S~ should be pos- 
itive and could be large enough to explain the NuTeV 
anomaly. 
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